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Abstract 


Given a linear time Invariant multivariable system = Ax^ t Bu_, ^ = Cx 
with m inputs and p outputs, Davison C63 has shown that p closed loop poles 

of the system can be preassigned arbitrarily using constant gain output feed- 

A A AAA 

back provided CA, B3 Is controllable. This paper shows that if CA, B, C3 is 
controllable and observable, and Rank B = m. Rank C = p, then max (m,p) poles 

of the system can be assigned arbitrarily using constant gain output feedback. 

Further, It Is shown that In some cases more than max (m,p) poles can be 

arbitrarily assigned. A least square design technique Is outlined to approximate 

the desired pole locations when It is not possible to place all the poles. 
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Introductlon 

The design of linear multivariable control systems using output feedback 
has attracted the attention of several authors. Cl “33. There are two ways of 
approaching this problem. The first method consists of estimating the states 
of the system using an observer and use these states In the subsequent design. 
In the second approach, either static or dynamic feedback of the output is 
used directly in the control problem and this view Is adopted here. 

Consider a linear time- invar i ant multivariable, system 
x_ = Ax^ + Bu_ 

( I ) 

^ = Cx_ 

where x is an n vector of states, u Is an m vector of Inputs and y_ is a vector 
of outputs. It is well-known that the problem of pole assignment using state 

A A 

feedback Is equivalent to the controllability of the pair (A,B) [43. Pole- 
shifting techniques for multivariable systems using static feedback has been 
studied by Retallak [53, Davison [63 and others. It has been shown by Davison 
that If (A, B) Is control lable, and If Rank C = p, then p poles of the system 
can be arbltrar I ly 'p laced using output feedback. This paper shows that given 

AAA ’A A 

(A, B, C) controllable and observable. Rank B = m and Rank C = p, then at 
least max (m,p) poles of the closed loop system can be arbitrarily placed using 
output feedback. 

The above result may be used In designing systems for high integrity in 
the event of failure of transducers and/or acuators [73* Due to failure there 
may be loss of inputs and/or outputs In the system. In such an event' this 
design takes advantage of the unequal number of Inputs and outputs so as to 
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assure no loss In pole assignability. 

Theorem 

A /\ 

Given the system (I) with Rank B = m <_ n and Rank C = p <_ n, then a 
linear feedback of the output u_ = Ky, where K is a (mxp) constant gain matrix, 
can always be found such that max(m,p) eigenvalues of the closed loop system 

can be be made to take preassigned (complex eigenvalues occurring in 
conjugate pairs) values. 

Proof : 

Let (A|,A 2 » — ,A n ) and-(p|,p 2 » — ,P n ) be the eigenvalues of the open- 
loop and closed- loop system respectively. 

We have 

open loop characteristic polynomial = |sI-A| = (s-A|)(s-A 2 ) (s-A^ 

(2) 

and closed loop characteristic polynomial = |sl-A+BKC| = (s-p ( )(s-p 2 ) (s-P n ) 

(3) 

Then 

| S 1-A+BKC | = de+ j-j. + g K 5 ( SI _ A ) -1 ]. (4) 

! SI-A | 

Choosing K = fd T whe#e f is a mxl (column) vector and d T is a Ixp (row) vector, 
and using the Identity det Cl + MN] = det Cl + NM], equation (4) becomes 

| — : A+BKC | = I + d T C (sl-A)~ 1 Bf 

T" _ I - _ I _ | - 

= I + d CT (sI-T AT) 1 T Bf (5) 

= I + d T C (SI-A)" 1 Bf 

1 - 1 * - 1 * 

where C = CT, A = T AT, B= T B and T is a nxn nonsingular matrix. 
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For clarity, the theorem Is Initially proved for the case of distinct 
eigenvalues of A and the multiple eigenvalues of A are considered In the latter 
half of the proof. 

Distinct Eigenvalues 

In this case equation (5) gives 


1 sI-A+BKC 1 
jsI-Aj 


n a , 


i-t (s_A i > 


(6) 


The value of oij depends on the closed loop eigenvalues (pj, — , p n >. 
From (5) and (6), 


T l n a. 

d'c (sI-A)' 1 Bf = l 

1=1 15 I 

Choosing T as a modal matrix equation (7) becomes 

t i n 

d C (sI-A)” Bf = l “i 

i=l (S-Xj) 


where A= diag. (X j ,A n ). 

Let c' be the i +h column of C and bj be the I th row of B. Then, 

a ( = d T c*b.f - ,1 = 1,2, — , n. 


(7) 


( 8 ) 


(9) 


Case ( i ) 

Let p>m i.e. more outputs than inputs. Choose fj such that bjf = 6j t 0 
i = 1,2, — , n. This can always be done since bj ^ 0, for controllability. 

Hence, d^c* = ot j/6 j 1 = 1, , n. (10) 

This gives c"*"d = a (ID 

where 


a_ = col Cci|/6|, c^/^* , 

Now, let Cp be the matrix made of the p independent rows of C^ and the 
corresponding subset of a. Then, 

d = C"' a 
P “P 


( 12 ) 
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Thus (dj, d 2 , — > d ) can be chosen corresponding +o the p desired pole 
locations. Once this is dene the remaining (n-p) poles are fixed automatically. 
Case ( i i ) 

Let m > p i.e. more inputs than outputs. 

Choose dj such that d^c 1 = y. ^ 0 1=1,2, , n. 

This can always be done since c^ I 0, for observability. 


Hence, b.f = Oj/yj 1=1, — , n. 

This gives Bf = a_ 

where a * col [a ( /6 ( , a 2 /y 2 , , ct n /yj. 


( 13) 
( 14) 


n 'n - 


Since the rank of B ism, there are m independent rows of B, B m> such that 

B f = a 
m -m 


where is the corresponding subset of ct. 

f » B -1 a 
m -m 


(15) 


Thus (fj, f^, , f m ) can be chosen corresponding to the m desired pole 

locations and the remaining (n-m) poles are located automatically. From case 
(i) and case (il) it Is evident that at least max (m,p) poles of the system 
can be assigned arbitrarily. 

Multiple Eigenvalues 

Let the eigenvalues of matrix A be A ( , > A w w ^ +h multiplicity 

n j, n 2> — , respectively. Choose T such that A = T AT has the Jordan 
canonical form with o blocks of respective sizes nj, n^, — , n^ and Aj, A^, — , 
A^ the correspond! ng eigenvalues. 


Now, we have 


sI-A+BKC I 
jsI-A 


(s-Pj)(s-p 2 ) — (s-p n ) 

rn n? n, , 

(s-A.) (s-A„) (s-A ) w 

I 2 0) 


(16) 
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CO 


where £ n. = n. 
i = l 1 

Equation (16) can be re-written as 


sfrA+grcl . | + J 1 _!l_ + ... + £ Jll 


n to 

a) a , 


sl-A 


( 18 ) 


i = l (s-A , ) 


1 = 1 (s-A ) 
a> 


The value of a'j ( 1 = 1, n^, j=l, io) depends on the closed loop poles 


(Pj» P 2 > "“P n )- 

From equations (5) and (.18), we get 


n i 


d C(sl-A) - Bf = l 


I 


1 = 1 (s-A,) 


+ l 


n co 

co a. 


(s-A )' 

CO 


( 19) 


(sl-A) has the quasi -d lagona I form diag D , , > J (i D where J, Is a 


n, x n, matrix of the form 


i 


i 


(s-A , ) 


(s-A j V 


( S-A f ) 


(s-A,) ! 


/ 1 xh,“l 

(s-A.) i 


(s-A,) (s-A,)' 


(s-A , ) 


( 20 ) 


Let C = [C 1 , C 2 , C w 3 and B = B 1 , Q Z , B^U 1 where cH is a pxn. matrix 

J 

and B 1 - 1 is a n^.xm matrix. Then it can be easily seen that 

d T C j J.B j f= - , 

J |£| (s-x >' 


,1 o2 


,co-iT 


j = l ,2, co 


( 21 ) 
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Fur + her It can be shown that 

orj = d T CC'jB'j + + + C J ' B J 3f 

"j n j 

<4 = d T CC , | + + C j ,B J ’ 3f 

J J 

j = l >2, a) (22) 

a J ’ = d T dB J ’ f. 
n J 1 "j 

Where cj Is the I^* 1 column of and B^ Is the 1^ row of B^* . 

In the matrix form equation (22) can be written as 
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Case (i ) 

Let p>m i.e. more outputs than Inputs and B^'f = 6^, i = I, — , nj and j=l, — a>. 
For controllability, every row of B corresponding to the last row of each 
Jordan block of A is linearly Independent [8] J.e., i 0, j=l,2, — -a>, are 

n j 

linearly Independent. 

Now, we can choose (f^fj, — ,f ) such that 

B J ' f = i 0, j = I , oj . (24) 

J J 

Substituting this In equation (23) we get 



II 12 2 

Where a = col Ca.,a„ ---a , a , - - - a , 

— 12 n . I n 0 
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Deflne a quasi -d iagona I matrix M, 


M 4 diagCM. ,M_, M 3 


where Is given by 


M .= 
J 


e i s i 


'j 


0 0 


0 ) 


S J 
n . 
J 


6 ^ 6 ^ 6 J ’ 0 


we have MC d= a 


(25) 


to i n i 

M is a nxn non-singular matrix since det M = gj (6 r ) t 0 by (24). Hence, 

T ' T 

Rank M=n and Rank MC = p. Let be the p independent rows of MC and let 

Op be the corresponding subset of a. This gives 


or 


C d = a 
P “P 

d = C"'a 

p -p 


(26) 


Equation (26) is similar to equation (12) and the rest of the proof follows 
as in the Case (i) of distinct eigenvalues. 

Case ( i i ) : m>p i.e. more inputs than outputs. 

Let d T C*j = «•{ 1=1,2, — , nj , j = 1,2, — u>. 

For observability, every column of C corresponding to the first column of 
each Jordan block of A Is linearly Independent [8j i.e. C'j ^ 0, j=l,2, — to. 

are linearly Independent. 

Now, we can choose (d^d^, — d^) such that d^C^ = 6^ / 0, j = 1,2, — to. 


( 27 ) 
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Subs+t+u+Ing this In equation (23) and defining a quas I -d iagona I matrix N, 


diagCN | , Nj, N j 


where hL is given by 


1 


! 0 


6^ — 5 J 

2 n . 

J 

i n 1 

J 


0 


I 


(t Is seen that 
NBf = a. 

w j n j 

N Is a nxn non-singular matrix since det N= j g [ > J f*0 by (27). 

Hence Rank NB=m and let B be the m Independent rows of NB and let a be the 

m — m 

corresponding subset of a. This gives 

B f = a 
m -m 


or f = B a . (28) 

m — m 

Equation (28) Is similar to equation (15) and the rest of the proof follows 
as in the case (II) of distinct eigenvalues. 

This completes the proof in the case of multiple eigenvalues. 

Nlature of The Design Equation (23) : 

In general, the output feedback gain matrix K = fd^ is obtained by 
solving the set of n non-linear simultaneous equations In (m+p) variables 

(d^d^, , dp , f j , f f m ) • However, in the proof" of the theorem either 

(d.,d„, — d ) or (f.,f„, f ) are selected arbitrarily and (23) is reduced 

■ to a set of linear equations and this assures at least max(m,p) poles 


can be placed arbitrarily. 
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In certain cases the non-linear nature of (23) can be exploited to assign more 
than max(m,p) poles of the closed loop system. 

Complex eigenvalues of the matrix A present an interesting situation. 

_ I * _i - * 

The Jordan canonical form A = T AT and the matrices B=T B and C = CT will 
then be complex matrices. However, K will be real since the complex columns 
of T - * and elements of o occur in conjugate pairs. 

In designing the control system using the theorem of this paper, max 

(m,p) poles of the system can be assigned as desired and this fixes the 

location of the remaining Cn-max (m,p)3 poles of the system. Let us call 

these poles the ’’dependent poles", B, where B = (p . ..., , p p ). 

' — fnaxirn^p/rj n~ I n 

In some cases, by taking advantage of the non-linear nature of (23) more than 
max(m,p) poles can be arbitrarily assigned and this reduces the number of 
dependent poles. However, nothing can be saida priori about the location of 
these dependent poles. Now, the would be a function of By minimizing 
a least square error criterion of the form 

^ T 1 9 

J = l q, [d C B f - o.AQ)T 
1=1 1 11 

subject to the constra 1 nts g ( £0 > 0 an approximate set of the desired closed 
loop poles can be realized. The weighting coefficients qj can be used to 
control the error between a pole In the desired set and Its corresponding pole 
in the approximate set. A somewhat similar approach has been suggested by 
Jameson C93 and Fallslde ClOU. 

If no satisfactory set of poles results from the least square error 
approach, e.g., closed loop poles unstable, then a dynamic compensator C3] 
would be necessary for pole-placement. It should be noted that some systems 
which need a dynamic compensator for pole placement using C33 can be made to 
attain any desired closed loop poles using only constant output feedback 



wi+h the above design method. 
Example I: 



The system Is controllable and observable with two unstable poles at I and 2. 
Also, m = 3 and p = 2. According to Davison two poles can be placed arbitrarily. 
However, according to the theorem stated In this paper three poles can be placed 
arbitrarily. Pearson would need a first order compensator to control all the 
poles. Here, it will be shown that by solving the equation (23) In Its non- 
linear form all the four poles of the system can be placed arbitrarily. We have 
d T C Bjf = a. I = I, 2, 4, 4. 

so, 


Vi ' “i 


f d = ot_ 

2 I 2 


Vz = “3 


<f l +f 2 +, 3 >d 2 “ V 


solving these equations with dj = I, we get 


d 


2 


V a 

a ,+a 


3 

2 



^2 = a 2 an< ^ ^3 


a^Ca | +a 2^ 
V a 3 
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(a 4 ~a 3 )/(a |+a 2 ) 

K =fd T = a 2 a 2 (a 4 -a-j)/(a| + a2^ 

a^(a i+c^/la^-a^) a 3 

with this choice of K all the closed loop poles can be placed at the desired 
locations. 

If the closed loop poles are desired at -l, -2, -3 and -5, then a^-7.2, 

0.2 - 14, a-j, = 0 and = 0.2. This gives f ( = -7.2, f 2 = 14, f^ = 0, d ( = I, 



This problem Illustrates the nature of (23) when A has complex open loop poles. 
The open loop poles are at I and - ^ ± j /3/2. If the modal matrix T and Its 
i nverse are chosen to be . 
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Theri A, B and C are given by 


A = T~’at = 


1 

0 

0 

1 ./3 

0 


| 

./3 

- 2 + J— 

0 

1 ,/3 

, B = T 'B= j 

~ 2 ~ 

J 2 

0 

0 

- 2 " J- 2 



J 2 _ 


and C = CT = 


1 I 

7 I . ./3 

2 2 + J— 


1 

2 - J 2 


If the closed loop poles are chosen to be at tj and I, then we have the set of 
equations 


j Cd ( + 2d 2 lf = 


(29) 


~ C-d , + d„3 - j Cd . t d«3 = a 


6 C " d , + d 2 ] + j 


2/3 1 

— — Cd. t d 7 3 = a 
2/3 1 2 3 


(30) 

(31) 


for the desired pole assignment. 


where a. =0, a„ = ■* 4- ^ — and a, = - ^ + -“J— 

I 2 2 2/3“ 3 2 2/3 

Equation (30) and (31) are complex conjugates and give the same set of equations 

in (d|,d 2 ,f). From equation (29) and the real and imaginary parts of (30), we 

get 


(d j t 2d 2 )f = 0 
(d j. - d 2 )f = 3 
(d | t d 2 )f= I . 


(32) 


Solving (32) with f = I, gives d| = 2, d 2 = -I and K = [2 -G. This choice of 
K gives the desired pole-placement. 
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Conclusions : 

It Is shown that for a controllable, observable system at least max (m,p) 
poles of the system can be arbitrarily assigned if Rank B=m and Rank C=p. Further, 
it is shown that in some cases more than max (m,p) poles can be arbitrarily 
assigned. A least square design technique is outlined to get an approximate 
set for the desired pole placement when It is not possible to place all the 
poles. The application of this technique for high Integrity design of Sikorsky 
SH-3D Sea King helicopter Cl I, 12] Is under study. 
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